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ABSTRACT 



The finite-amplitude standing wave behavior of air within 

several tubes was experimentally Investigated near standard 

conditions of temperature and pressure for frequencies below 3000 

cps. The sound pressure levels obtained were large enough to 

generate shock waves of 0. 10 atmospheres. The air in the tube 

was driven by means of a vibrating piston. The motion of the 

piston was measured by an accelerometer while the pressure at the 

rigid end was measured with a condenser microphone. It was found 

that the wave forms for different tubes would be the same if the 

following quantities were made equal: 1) the ratio of the driver 

acceleration to the acoustic attenuation constant, and 2) the phase 

angle between the acceleration and the pressure at infinitesimal 

amplitudes. Values of the attenuation constant were determined by 

several different methods including measuring the decay of 

pressure after clamping the driver piston and by determining the 

ratio of acceleration to pressure. Observed attenuation constants 

-4 —1 

were of the order of 10 cm . The resonant frequency for 
infinitesimal amplitude was observed to be the frequency maximizing 
both the average rectified pressure and the shock strength only for 
the stronger shocks. However, for weak shocks the shock strength 
became a maximum at lower and lower frequencies as the pressure 
decreased. This tendency, as Indicated by the amount of second 
harmonic distortion, continued in the weak finite amplitude region. 
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1. Intrcxluction. 



The purpose of this research is to investigate the parameters 
that effect the behavior of high intensity standing waves in air 
contained in tubes. In many practical instances the sound 
pressure levels within cavities are found to be of such large 
magnitude that infinitesimal amplitude theory no longer applies. 

The investigations reported here were carried out in air-filled 
tubes each terminated with a rigid boundary at one end and driven 
by a vibrating piston at the other. 

The acoustic pressure at the fixed end (x = L where L is the 
tube length) predicated by the theory of infinitesimal amplitude 
waves with no dissipation is given by the equation; 



The resonant frequencies are given by kL = nTT, and at resonance 
the predicated pressures are infinite. 

If dissipative terms are introduced into the acoustic wave 
equation, then the pressure at x = L becomes: 



Where 

A = acceleration amplitude of 



H S /‘n KL 



the driver piston at x = 0 
k = wave number 6>/c 
c = local speed of sound 
Po = density of the gas 



^ A 
K 




1 




4 







Where o(, is the attenuation constant depending on viscosity, heat 
conduction, molecular relaxation of the gas, and the wall effects. 
This solution is still not complete: although it eliminates the 
singularities at resonance, it is still a linear theory predicting 
continuous pressure-time behavior near resonance. This Is an 
inadequate description since experimenters have observed that 
shock waves are present in the tube near resonant frequencies. 

In the region between infinitesimal amplitude waves and 
shock waves seyeral things are happening. As the amplitude 
increases the sinusoidal pressure wave begins to distort as a 
result of increased inportance of non-linear effects in the gas. 
These non-linearities generate harmonics of the fundamental 
which lead to an apparent distortion of the observed waveform. As 
the strength of the fundamental is Increased, the distortion grows 
until a shock wave finally developes. 

Most investigators of finite amplitude effects have confined 
their studies to traveling waves (2, 5, 8, 19, 20). There have, 
however, been a few investigators of finite amplitude standing 
waves (1, 4, 17), 

In this study a tube terminated at one end by a rigid piston, 
and driven at the other by an oscillating piston was used to 
investigate finite amplitude standing waves. The use of a piston 
as a sound source in a tube is no new concept (11, 12, 18). How- 
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ever, most pistons in the past have been driven by some sort of 
reciprocating mechanical driver. The space program has created 
the need for vibration exciters to test missile parts, and out of 
this need has come the machine used in this experiment as a 
piston driver. Reichwein (13) first used a driver of this type with 
a capability of 16 g acceleration to drive the piston in a tube. 

The tube was terminated by a plug containing a condenser 
microphone to monitor the pressure at the terminal end. He 
investigated threshold behavior of periodic shock waves in 
resonating gas columns and reached the conclusion that in a 
given tube the waveform was a function only of the acceleration 
of the driver piston and the "distance" off resonance as measured 
by the phase angle between the shock front and the piston 
acceleration. 

Sanders^ later observed that the waveform depended markedly 
on the "tightness" of the tube enclosure. He found that if there 
were leaks at either end of the tube shock waves would not 
develop until the acceleration was much higher than for a "tight" 
tube, and that the waveform in a rectangular cross section tube 
was different from that in a round cross section tube. It was 

^Sanders, J. V. , private communication. 
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postulated that the tube attenuation was an additional parameter 
affecting shock wave development. 

The attenuation of infinitesimal sound in a tube has been 
treated by a number of people both theoretically and experi- 
mentally (6, 7, 13, 16, 21). It turns out that the attenuation 
constant in a tube is proportional to the square root of frequency 
due to the wall effects, plus other terms attributable to free 
space attenuation, etc. In this experiment tube effects are 
about four orders of magnitude greater than the other effects. All 
of the investigators have been able to verify theory to reasonable 
accuracy, but all have measured an attenuation constant higher 
than the theoretical one. Different explanations have been 
advanced for the discrepancies, but none have been proven. 

The experimental set-up used in the present experiment is 
essentially the same as that used by Reichwein. However, 
the 16 g vibration exciter was replaced by one of 50 g capability. 
The original experimental set-up was further modified by the 
addition of a water jacket on the tube to stabilize tube temper- 
ature. This arrangement enabled the experimenters to study 
standing waves from the infinitesimal amplitude region to the 
finite amplitude shock region without any equipment changes. 

2. Equipment Description 

Oscillator The Hewlett-Packard Model 205 AG Audio Signal 
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Generator was found to be a very stable oscillator. In the 
frequency range of 20-20,000 cps there is a drift of two per cent 
or less, and if allowed to warm up for thirty minutes the drift was 
less than 0.01 cycle in the frequency range of interest. 

Piston driver The MB Model CIO VB Vibration Exciter was 
selected for the wide range of accelerations attainable. The 
electrically driven piston is preferable to a mechanical or electro- 
mechanical device because the acceleration of the electrical 
driver can easily be controlled and there is enough power to 
drive a 0.5 lb load with 50 g acceleration or a 30 lb load with 25 
g acceleration. The dynamic range for the loads imposed was 
from a maximum of 50 g acceleration to a minimum of 0.5 g where 
60 cps noise made the signal to noise ratio unfavorable. 

Microphone The Brifel and Kj^^r Condenser Microphone 
Type 4136 was used for the pressure pickup. Since the 
microphone was embedded in the rigid end of the tube it was 
necessary to use a microphone of small diameter so as to 
produce as small perturbation of the boundary conditions as 
possible. For this reason the 0.25 in microphone was selected. 

I 

In addition, for the investigation of finite amplitude waves a 
microphone was needed with a relatively short rise time. The 
rise time of this microphone is 4,5 micro seconds and the 
frequency response is essentially flat from 25 cps to 80,000 cps. 
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Freouencv counter Hewlett-Packard Electronic Counter 



Model 522B was used for frequency determination. The time 
interval unit on a ten period average was used for frequencies 
below 1000 cps and the frequency unit on a ten second count was 
used for frequencies above 1000 cps so that a 0.1 cps accuracy 
in frequency determination could be obtained. 

Accelerometer Endevco Model 2215 Accelerometer was 
bolted to the driver face plate. This accelometer had been cross 
calibrated by Sanders yielding a sensitivity of 4.82 rms mv/peak g. 

Oscilloscope A Tektronix Type 565 Dual Beam-Dual Time 
Base Oscilloscope was used with plug in units 2A60 and 3A72. 

With these plug in units three tracers, could be viewed simult- 
aneously. The dual time base permitted the expansion of any 
portion of the waveform, and in particular the shock front for 
detailed inspection. The delay feature incorporated in this 
scope also allowed phase between two traces to be read directly 
from a dial. 

Tube enclosure Several tube arrangements were used to 
investigate a set of different boundary conditions. The basic 
arrangement is shown in Fig. 1-1 and is considered to be the 
most practical set-up since the tube length can be readily varied. 
The tube was of aluminium, length 5 ft. 11 in. , wall thickness 
of 0.20 in. , and an inside diameter of 1.75 in. Early in the 
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experiment it was found that driving the piston at large 
amplitudes caused the temperature of the air in the tube to 
increase thereby changing the resonance frequency. Therefore 
the tube was fitted with a double wall so that water could be 
circulated between the walls to stabilize the temperature. This 
water-jacket kept the temperature in the tube constant at 14.5°C 
by absorbing the heat from the inside of the tube and also 
shielding the inside from the heat produced by the electronic 
equipment in the room. 

The microphone was mounted in a movable piston eight 
inches long with its head flush with the face of the piston. The 
microphone piston and the driver piston were placed in the tube 
with "O" rings and vacuum grease to produce an air tight seal. 
The driver piston was connected to the face of the driver plate by 
means of a short threaded shaft. 
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3. 1 Theory of Infinitesimal Amplitude Plane Waves in a 

• m 

Dissipative Fluid 

The standard form of the linearized plane wave equation 
using particle displacement as the acoustic variable is 
^ j ^ A solution of this eauation is 

where 0 is the complex displacement amplitude of a plane wave 
of frequency to and wave length constant K, traveling in the 
positive X direction with velocity c; and D is the amplitude of a 
similar wave traveling in the negative x direction. The 
subscript tilde will be used to denote complex quantities. 

This solution and the wave equation for which it is a 
solution are based on the assumption that there are no 
acoustical energy losses, either in the medium or at the 
boundaries of the medium. Appropriate mechanisms have been 
developed to deal with these losses, some of which will be 
discussed in section 3.2. Assuming that these losses exist and 
that they are a function of the distance traveled by the plane 

t ' 

wave 'through the medium, it is possible to modify the plane 
wave equation to account for these losses. If the expression 
K = Ci>/c is modified to /^ - ^/c where C is a complex velocity 
and K is considered a complex wave length constant with real 
part K and Imaginary part -j then it can be shown that the 
plane wave equation may be re-written as ^ i C ^ ^ 
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with a solution = 3 ^ q ^ 

which may also be written ^ ~ 

This solution is similar to that for the unattenuated wave except 
that the displacement amplitude is attenuated in accordance with 
the term . The acceleration of the particles for infinite- 

simal waves may be found by taking the second derivative of 
particle displacement with respect to time, ^ 

The pressure of the wave can be expressed as 



= -yO ^ 

/w / O -V ^ 

^ X. 



= ^ 

(H -oc<y 



(j K icO Sc M)pe"V 



= [3 



e. e. — jj e ^ 



Let us now consider a plane wave of this form being 
propagated in a tube closed at one end and driven at the other by 
a piston moving with an acceleration Ae^^^ where A is the real 
acceleration amplitude. We take the piston rest position to be 
at X = 0 and the rigid end to be at x = L. 

Two boundary conditions are imposed: that the particles be 
stationary at x = L, and that they be moving with acceleration 
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at X = 0. 



At X = 0 
At X = L 



/I 






] 



C r ^ W L r% 

^^[Ber^e ^ J<^ 






Solving for B and D in terms of A 






r> — A 

kr ~ a>'‘^ re.****‘£*'^ 

When these values of g and g are substituted in the pressure 
equation, pressure becomes 

L ^«4 M-i. _ . ^ 



Evaluating this expression at x = L 



^ pg e-'^eju_i _. j- 






i d>^ 



[ilhh c^w cos Ki. ti CO iUt i'« 
The afriplitude of the pressure at the rigid end, p — 



P = 





/t 1*1 r j^s/nh^^L cos^ffL i-Ccs/i^o(L Sih^iriJ 

For the range of experimental values investigated eL is of 

-4 -2 2 
order 10 , K of order 10 , and L of order 10 . Since KL2J 1 

and olL«KL, the approximations can be made that sinh o(.L= o(L, 

cosh o(L = 1. Using this and the fact that << / ^ 



P 3 jL- T—J— 

K L (o(U C.CS i- SifO* KLj 
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^ 0 






f, 








I 



m 





This equation for P is correct to first order in o<v L. Finally 
may be expressed as where 0 is 



-I 






— tern * ^coih o(L ~tan 



To investigate the acoustic behavior of the system as the 
frequency or length is varied, the derivative of P with respect 
to either variable is taken and set equal to zero. 

For P = j^(o(L cos ^ 

setting k£. -.0 leads to 

,.f<C0S fU [ -K(e(L)^S>'n KL L CoS KL -O 

When cos KL = 0, P is a minimum. The requirement for cos KL = 0 

is that or L = ^ 

K 



Pm I n " ^ 

K 



o r P* ^ ^ 

Qrt - TT 

This leads to the conclusion that the same value of minimum 



pressure will be observed at lengths of the tube corresponding 
to '"k) >> where m is an integer greater tfian or equal 

/r * 

to one. The maximum value of pressure /occurs when-' 

,that is ~ n TT to first order in ot L. 

wD-K'-i'] 

The value of P is A 

-^rnr 

Taking the partial derivative of P with respect to frequency 
and setting it equal to zero leads to a determination of the 
resonances and antiresonances of the system, but the result is 
more complex since the parameter is also a function of 
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frequency but not of length. 

Fig. 3-1 illustrates the dependence of P on length of the 
tube at a fixed frequency, and Fig. 3-2 illustrates the 
dependence of P on frequency at a fixed length. The frequency 
dependence of << is taken into consideration in Fig. 3.2. 

3.2 Attenuation of Sound Waves in Tubes 

Acoustic waves being transmitted in a medium must undergo 

losses of acoustic energy. These losses are of two types, 

losses due to dissipation in the medium itself and losses 

associated with conditions at the boundary of the medium. 

Losses in the medium may be viscous losses, heat losses and 

losses due to molecular exchange of energy. For dry air at 

standard conditions these losses in the medium may be 

represented by an attenuation coefficient of the approximate 
“13 2 

value 2x10 f where f is the frequency of the acoustic wave 
in cycles per second and the units of «<. are cm”^. For the size 
of tube and frequency range dealt with in this experiment these 
losses in the medium are several orders of magnitude smaller 
than the losses due to boundary conditions and may be safely 
ignored . 

Losses associated with boundary conditions may arise from 
the direct absorption of acoustic energy by the walls. The 
requirement of zero particle velocity at the walls leads to a 
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velocity gradient from the center to the walls and thus to 



dissipative forces which are related to the shear viscosity of 
the medium. In addition, since temperature is essentially 
constant at the walls, there is a tendency for heat to be 
conducted to the walls from the medium during a condensation 
and from walls to medium during rarefaction, resulting in a 
dissipation of acoustic energy. 

The theory of attenuation in a tube of circular cross section 
has been developed by Helmholtz, Klrchoff and Lord Rayleigh [l^. 
A simplified exposition of this theory may be found in 
"Fundamentals of Acoustics" by Klnsler and Frey The 
attenuation in a tube due to boundary effects is shown to be 

oL = -k-r-7^ ^ - where a is the radius of the tube, ^ is 
^ ' SI 

an effective coefficient of shear viscosity, to' is radian 

* 



frequency, c is free field speed of propagation in the medium, 
and ^ is density of the medium. The speed of propagation in a 
tube is slower than in an infinite medium and is given by the 



effective shear viscosity rather than an actual coefficient of 
shear viscosity is a technique developed by Kirchoff and later 
Rayleigh to take into account the effect of heat conduction to 
the walls. The relationship between^ and ^ is 



expression 




The use of an 




where k is the thermal 
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conductivity of the gas, y is the ratio of specific heats, and Cp 
is the specific heat at constant pressure. Investigators have 
stated that in a tube of rectangular cross section the attenuation 
due to wall effects is directly proportional to the ratio of the 
perimeter to the cross sectional area , in place of the 
reciprocal of the radius in the case of the cylindrical tube |^, 16 j. 

The derivation of the thpory for attenuation in a rectangular 
tube, is complicated by the fact that it must be accomplished 
using rectangular coordinates, and an extra dimension must be 
considered which \was not necessary in the cylindrical tube due 
to its circular symmetry. The simplest means of approach to 
the rectangular case is to consider the case of a plane wave 
being propagated between two parallel flat plates in a direction 
X parallel to the plates. Consider an element of volume 
between these plates of length dx in the direction x, of unit 
width in the direction parallel to the plates and perpendicular 
to the direction of propagation, and thickness dy in the 
direction perpendicular to the plates. This element, of volume 
dxdy, has exerted on it a force due to the gradient of pressure 
in the X direction of magnitude . In addition, due 

to the boundary condition of zero particle velocity at the walls 
there will be a particle velocity gradient in the y direction. 

Since the viscous force on a surface of area dx is ^ m 
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where ^ is the coefficient of shear viscosity, the net force due 
to the viscous forces on the element of volume will be 






If these forces are equated to the 
mass of the medium contained in the volume element multiplied 
by the acceleration , a partial differential equation in JUL is 
obtained as follows 



~l£- ^ + a 

If the motion is assumed to be sinusoidal and of frequency u> 
then yb may be expressed as and ,44 as 

<j so that ^ may be replaced by Jo; , 

. X. “I 

u 



Then 


K 


= [o> 


CJ - 






or - ' 




- r 


, 1 ■s ■ 


"I 




yy 






^-K 


J y 




where 




^ j /o g> 


and 




K = ' 






A solution of the equation is ^ 

Applying the boundary condition that L/* = 0 at Q/2 and -CJ/2 



+-/A 



where a is the distance between the plates 

A _ P~ 



u 









The average value of \J over the distance a, since U is 
symmetrical about y = 0 may be found by integration as follows 

15 
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Since K = 







- j >f4 

X. 








Rjr » /o which is the case for the experimental 

! 9,M 

situation , 4< I and can be neglected. 

(7 then is approximately Q. 

0 - 



Solving for 

which is approximately 






Since 



/r 



X . 



K ^( 1 - 0 /^ 



1 (i+j)l y 

From this point the derivation follows exactly that of the 
cylindrical tube. The reactance term j^6>is the normal 
reactance term associated with and the effect of 

viscosity is to Introduce an additional reactance term ^ 

so that the effective density of the viscous medium is a ^ 



The effect of this effective density on the medium leads to the 
expressions for c’ and ^ previously given for the cylindrical tube. 
The conclusion is that the attenuation of a plane wave between 
two parallel plates is Inversely proportional to the distance 
between the plates, Just as it is Inversely proportional to the 
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radius in the cylindrical tube. Generalizing to the case of a 



rectangular tube, it is apparent that the attenuation due to the 
two sets of parallel plates making up the walls should be 
additive if comer corrections are negligible ^ The attenuation 



perimeter to the cross sectional area. This relationship agrees 
with the statement previously made that attenuation is directly 
proportional to the ratio of perimeter to area. (Fig, 3-10) shows 
the comparison of obtained using this theoretical development 
with the vK actually measured in the rectangular tube.) 

In the special case of standing waves in a cylindrical tube 
Parker [l3] states that the attenuation due to the wall effects 
should be increased to take into consideration the losses at a 
rigid end. His relationship is o(.t -o(wv where R is 

the tube radius and L is the tube length. This additional 
attenuation is of the order of three per cent for the tube size 
used in the experiment. 

3.3 Investigation of Attenuation in Tubes 

The investigation of attenuation in the infinitesimal 
amplitude region carried out in this study consisted of 
measuring attenuation by five methods suggested by the theory 



should then be proportional to where a and be are 

the dimensions of the cross section of the tube, /jr ■♦"•4^ 



be expressed as 




which is one half the ratio of the 
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previously developed. These methods were then checked for 
consistency and finally compared with the attenuation predicted 
by theory. 

The first method of measuring attenjiation Involved measuring 
phase changes of pressure with respect to acceleration by 
varying frequency in the vicinity of resonance. In^ the expression 
^ C ' where 

^ n 1 -han HL ^ 

0 represents the phase angle between the acceleration of the 

piston and the pressure at the rigid end of the tube. If the 

( 

same simplifying assumptions made previously are applied to 
the equation for 9 it may be solved for o< with the result 

^ s — where f is the frequency for which 0 is ^ 

L o 2 

and Zif is the change in frequency from f^ necessary to 
obtain some phase angle other than ^ . The method becomes 
more accurate as larger values of ^f are obtained, since A f is 
a small quantity of the order of a few cycles per second. The 
experimental procedure used was to vary the frequency about 
the center frequency f^, read values of 0 on a dual beam 
oscilloscope, calculate values of attenuation for each set of 
readings and compare consistency. The accelerometer and 
microphone voltages were displayed as vertical inputs to a 
dual beam, time shared CRO. The Lissajous technique was 
also tried but discarded since it introduced additional labor 



18 



into the method without increasing accuracy. Phase meters 
were also tried but the relatively low Input impedances of 
these meters significantly loaded the circuits and introduced 
phase changes. Throughout the experiment an attempt was 
made to keep electronic circuitry to a minimum of complexity, 
since amplifiers and filters introduced phase changes which 
were both amplitude and frequency dependent in an unpredict- 
able manner. In the early stages of the experiment it was 
discovered that the value of f^ was extremely sensitive to 
temperature changes in the medium, a one degree centigrade 
change corresponding to a 0.34 cyle per second change in Iq 
at 200 cps. Since Af was of the order of one cps, this was a 
signiflpant source of error. The problem was eliminated by 
water jacketing the tube to stabilize the temperature, after 
which drift in the value of fQ became negligible. Fig. 3-3. 

The expression r from which — 

provides a simple and rapid means of determining . Since KL 
is apprc»dmately nTl’the parameters to be measured are accelerf* 
ation and pressure, ^9 being obtained from a handbook as a 
function of temperature. Assuming proper calibration of both 
accelerometer and microphone, this method should produce 
accurate results. Since KL may be made any one of many 
values of n within the frequency range of the equipment. 
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values of «< for many different frequencies may be obtained with- 
out varying the length of the tube. Fig. 3-4. 

The e< obtained by measuring P and A should be the same at 

all the modes n obtained by varying length and holding frequency 

fixed. The product KL P ^ should be a constant for all n for a 
^ max 

given A. The experimental results of plotting 3is ordinate 

against n for a fixed value of A was not the expected strai^^ht line 
but a curve which asymptotically approached a horizontal line for 
high values of n. Fig. 3-7. The explanation of this phenomenon 
comes from a consideration of the losses making up o( . The 
theoretical value of o( was derived assuming a boundary condition 
of a rigid end at x = L. If there are leaks around the gaskets 
there will be significant increases in the value of attenuation 
measured. By using a single parameter o< the losses at the ends 
have been treated as though they were distributed over the whole 
length of the tube . A more appropriate treatment would be to treat 
attenuation of pressure amplitude for a wave traveling down a tube 
and striking an end as attenuation of the form 
where ^ is the loss in amplitude due to the end. If the 

length of the tube is increased the attenuation due to the end does 
not increase but the attenuation due to the walls is increased. As 
tube length becomes greater the losses at the end become less 
significant. Thus at higher values of n the end losses will be 
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proportionally less significant and the measured values of •< 
should approach asymptotically the value of «< for ideal end 
conditions. Fig. 3-8. 

By comparing the value of o< for n = 1 with the value of 
for a higher value of n at the same frequency an indication of the 
tightness of the system is obtained. As a result of these 
comparisons a greatly improved rigid end piston was designed and 
installed. With this piston the system became so leakprobf that it 
was difficult to change the tube length. It took several minutes 
in each case for pressure equilibrium to be established after a 
change in tube length. 



Pmin Is a function of frequency but not of length, that is, for a 

given value of frequency the values of P . obtained by varying L 

min 

through different values\df KL = (m-l/2)IT' will all be the same. 
This conclusion was verified within the limits of experimental 
accuracy. 



If the expressions for P_,„ and P . are combined, the result 



which Pjjiax obtained. Since o( may be found using a ratio of 
pressures this method of determining is considerably less 
sensitive to calibration errors. One objection to the method is that 
at the small values of ^min involved, system noise becomes significant. 



The expression *■ ^ leads to the conclusion that 

K 




is obtained, where L is the length of the tube at 
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An additional method of determining is to determine the 
frequency or length change necessary to reduce the value of P to 
0. 707 of its maximum value, that is to reduce it be three decibels. 
In this "bandwidth" method 

/pioox ^ -7- C0S>fL')^ ^ 

In the vicinity of resonance KL may be assumed to be some small 

value away from nTT and «>< may be assumed constant 

2 

for the small range of frequencies involved. Sin KL may be 

approximated by or, with L held constant, by 

2 

Cos KL may be approximated by 1. Then 

^ f (Ak7 

o< A K 

In practice it will be simpler to measure the change in frequency 
between down 3 db points on both sides of the resonant peak. In 

terms of frequency r TT A-f where ^ f is defined as the 

c 

frequency change necessary to go from one of the half-power points 
to the other, and c is the speed of propagation in the tube. This 
relationship was used by Parker [l3] to determine attenuation in a 
tube filled with oxygen. 

The bandwl,dth method is a simple method with the advantage 
of employing a ratio of pressures, thus eliminating calibration 
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errors. Since the f in the above expression is only of the order 
of a few cycles per second, and the determination of the driver 
piston frequency was only to the nearest tenth of a cycle per 
second, this method was not particularly consistent. A more 
accurately controllable system and a more accurate frequency 
counter would improve this method considerably. Obtaining the 
half-power points accurately on a voltmeter is also difficult and 
not consistent. 

A method of measuring o< independent of the previous methods 
is to measure the attenuation of the plane wave in the tube after 
the driving piston is stopped. When the power to the piston is cut 
off the piston stops in one or two cycles , leaving the plane wave 
to decay in a tube with rigid stationary ends at x = 0 and x = L. A 
picture of this exponential decay can be taken with a Polaroid 
camera mounted on the oscilloscope with the sweep time base 
adjusted to display 3 or 4 time constants of the decay. (Figure 3-5). 
If the magnitude of the pressure is plotted againipt time on semi- 
log paper a straight line is obtained, from which a decay constant 
is calculated. (Figure 3-9). This decay constant is in units of 
sec ^ but it may be converted to c{ in units of cm ^ by dividing 
it by the speed of propagation in the tube. 

The objections to this method are first that it is only 
appropriate for KL = 1 fT , since for higher values of n beats occur 
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which make analysis of the exponential decay more difficult. 

(Figure 3-5, 3-6). Second, it is slower. Third, non-linearities 
in the display system of the oscilloscope become significant. 
Fourth, the different end condition at x = 0 may result in a slightly 
different value of «<, than existed when the piston was moving. 
Even with all these objections, this method proved quite 
consistent and on a par with the other methods. 

The investigation of attenuation in the acoustic region in a 
tube of rectangular cross section was carried out in a similar 
fashion as in the round tube. The equipment used made it 
impossible to vary the length of the tube more than a few centi- 
meters. All methods of determining o( could be used with the 
exception of the method involving the ratio results 

of the four methods of obtaining attenuation used on the rectangular 
tube were consistent and in agreement with theory. Figure 3-J.O. 
Table 1 illustrates the comparison of theoretical values of 
attenuation for both the cylindrical and rectangular tubes with the 
attenuation measured by the five methods discussed. These 
results illustrate the consistency that can be obtained by making 
measurements under identical conditions. 
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4. Investigation of Finite Amplitude Waves in Tubes 

As discussed in the introduction, when a piston is 
sinusoidally driven with large amplitude in a rigid ended tube, 
the pressure at the end of the tube is not sinusoidal when the 
frequency is near resonance. It is observed that the pressure 
wave at the rigid end of the tube contains a discontiniulty . A 
pressure discontinuity of the type observed is characteristic of 
a shock wave. 

Saenger and Hudson (1960) developed a finite amplitude wave 
theory based on the assumption that the solution to the problem 
is the sum of a continuous part, which is calculated as a power 

series, and a discontinuous part, which satisfies the ahock 

/ 

conditions. This theory is that used by Reichwein (1962| in his 
Investigation of threshold behavior of shock waves in a tube using 
similar equipment to ours. Saenger and Hudson also did some 
experimental work and reported fair agreement with their solution 
based on a minimal number of observations. Two other 
theoreticians in the field of finite amplitude waves are Betchov 
(1958) who also constructed his solution from an assumed sum 
of a continuous part, and a discontinuous part, and Chester (1962) 
who refines a collection of the previous theories. 
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Effects of attenuation on waveform Figure 4-1 is a montage 



of oscillograms illustrating the growth of shock waves as a 
function of piston acceleration, and phase angle away from 
resonance in two tubes of different cross sections. The resonance 
frequency in both tubes was 89 cps, and the phase angle between 
pressure at the rigid end and acceleration of the piston was 
measured in the infinitesimal amplitude region with pressure 

having the leading phase. The normalizing factor was 81.3 

2 2 
m/sec in the round cross section tube and 93.8 m/sec in the 

rectangular cross section tube, and was the acceleration necessary 

to produce an undistorted wave of 1.57 x 10"3 bars rms pressure. 

It is clear from these oscillograms that there is an intermediate 
region between the infinitesimal amplitude and shock regions. 

This intermediate region is characterized by an increasing 
distortion of the infinitesimal amplitude wave as the acceleration 
of the driver is increased until a discontinunity appears and will 
be called the weak finite amplitude region. Saenger and Hudson 
predict a sudden appearance of the shock wave which was not 
observed. Instead there is such a gradual growth of distortion 
that there was considerable disagreement between observers as to 
when a shock was actually present or to when the wave was 

I 

merely badly distorted. From the rise times of the measuring 
instruments it was calculated that the rise time of the fully 
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developed shock wave was less than 4 micro seconds. Saenger 
and Hudson further predicted that a shock could not exist below a 
piston acceleration which is ten times larger than that observed 
at 100 cps. 

The fact that it was necessary to water-jacket the tube in 
order to keep the resonant frequency the same for high amplitude 
work yielded some experimental evidence that the air temperature 
in the tube was rising as predicted by Saenger and Hudson. The 
tube resonance conditions were not observed to change when 
infinitesimal ampljLtude waves were present in the tube. 

With the available frequency counter, frequency could not be 
determined with an accuracy greater than 0. 1 cps. This presented 
difficulties when it was required to reproduce measurements at a 
given frequency. Phase angle, wh^ich can be measured to 1.0°, 
is a more sensitive p>arameter than frequency as can be seen on 
the montage where two cps correspond to approximately 60° phase 
angle near resonance. Phase angles were adjusted by keeping the 
tube length constant and changing the frequency of the driver 
until the desired phase angle was observed. It was necessary, 
however, to measure phase angles in the infinitesimal amplitude 
region because, if the upgolng crossing of the axis is chosen as 
the point at which to measure phase, this point is not fixed as the 
amplitude of the piston is increased to the point of pressure wave 
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distortion; i.e. , the weak finite amplitude wave phase angle must 
be defined in a different manner than the infinitesimal amplitude 
phase angle and will not be defined in this report. The definition 
of phase angle used in the shock region is the phase angle 
between the shock front and the upgoing axis crossing of the 
piston acceleration. 

Data for the montage were limited to 60° off resonance by the 
acceleration capability of our equipment since observations that 
were made at phase angles further from resonance revealed no 
shock wave development at the maximum acceleration of the 
equipment. If more powerful equipment had been available shock 
waves might have been observed at frequencies further from 
resonance, but just how far off resonance shock waves will 
develop cannot be reckoned. From the shape of the montage it 
might be expected that if unlimited acceleration was available 
shock waves could be produced for any frequency away from 
resonance; however, the above mentioned theories predict the 
shock condition only at resonance and frequencies near resonance. 

The similarity of the waveforms in the two tubes is remarkable 
since the tube cross sections were so different. The reason for 
this similarity lies in the process olj normalizing the acceleration. 
In the Infinitesimal amplitude region pressure is proportional to 
the acceleration of the driver and inversely proportional to the 
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attenuation constant. (Refer to section 3.3) If the acceleration 
required to produce a given infinitesimal amplitude pressure is 
taken as anormalizing factor in one tube, and the acceleration 
required to produce the same pressure in a second tube is taken as 
the normalizing factor in that tube, then the effect of the different 
attenuations in the two tubes is negated by this normalization 
process. Even though the attenuation measurements made in the 
infinitesimal amplitude region do not apply in the finite amplitude 
region it is remarkable that the wave forms in the different tubes 
are similar in both regions when the accelerations are normalized 
in the infinitesimal amplitude region. Had the acceleration ipelf 
been used as a parameter (as predicted by Reichwein) the wave- 
forms would not have been in such close agreement. The 
attenuations within the tubes used were nearly the same, but 
tubes with quite different attenuations can be expected to produce 
identical waveforms if the accelerations are thus normalized. 

Weak finite amplitude region It was observed that visual 
distortion of the pressure wave appeared when the magnitude of 
the second harmonic was 10 per Cjsnt that of the fundamental. A 
wave analyzer was used to measure the second harmonic as the 
acceleration of the piston was increased until the magnitude of 
the second harmonic was 10 per cent that of the fundamental. 
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These measurements were made in a round cross section tube and 



in a rectangular cross section tube and the results are plotted in 
Fig. 4-2 and 4-3. The same procedure was repeated for 20 per 
cent second harmonic and for the threshold of shock. No 
quantative data were plotted for the threshold of shock, however, 
since the point at which the shock appears fully developed depends 
on the opinion of the observer. It appeared as though the threshold 
of shpok was near the point of 30 per cent second harmonic, but 

•I 

data taken for this region are not presented since the definition of 
shock becomes more obscure as frequency is taken away from 
resonance. The wave analyzer used introduced a phase change 
which was a function of frequency and amplitude which precluded 
determination of the phase relationship of the second harmonic to 
the fundamental. It is interesting to note that for a given 
acceleration or pressure the maximum distortion occurs at a 
noticably lower frequency than the infinitesimal amplitude 
resonant frequency. 

Shock region Reichwein investigated the parameter 
defined by Saenger and Hudson to be the phase angle between the 
maximum displacement of the piston out of the tube and the 
arrival of the shock wave at the piston. He measured velocity of 
the piston whereas we monitored piston acceleration resulting in 
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our phase angle parameter being 90® out of phase wif,h his ^ 
parameter. Reichwein drove his system at maximum amplitude, 
set his desired ^max by varying frequency or length of the tube, 
and then reduced the amplitude down by steps until the amplitudes 
of the piston acqeleration became so low that shocks could no 
longer be observed. The value of ^max for which no change 
in A occurred as amplitude decreased was not Amax = 0 as 
might be expected but was some value in the vicinity of A max = 
-10*^ dependent on the resonant frequency. He further predicted 
that if amplitudes of acceleration were increased above his 
maximum available acceleration, that A would remain constant. 

We repeated the same procedure using a maximum acceleration 
which was approximately three times that of Reichwein. The 
curves obtained, and presented as Figure 4-4, have the same 
shape as those presented by Reichwein. The ordinate of this 
graph is not A., but the phase angle between the shock at the 
rigid end and the upward going axis crossing of the piston 
acceleration. This phase angle is equal to A+ 90°. For 

O 

accelerations up to 200 m/sec (the maximum used by Reichwein) 
our results were in agreement with his, and as the acceleration 
was Increased further the predicted trend was, followed until about 
400 m/sec^. For accelerations above 400 m/sec^, Reichwein' s 
prediction that phase angle would remain constant could not be 
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verified due to piston acceleration distortion problems. 



Resonance in the finite amplitude region The definition of 
resonance in the infinitesimal amplitude region is that frequency 
at which the maximum rms pressure is obtained for a constant 
piston acceleration, and it has been shown that this occurs when 
the phase angle between pressure and acceleration is 90®. At 
the frequency of infinitesimal amplitude resonance, for all but 
weak shocks, the phase angle between the shock and the 
acceleration is 90®. If the piston acceleration is held constant 
and the frequency changed, it was observed that both the shock 
strength and the average rectified voltage (as measured on a VTVM) 
decreased. This is in complete harmony with the definition of 
resonance in the infinitesimal amplitude region. If the shock 
strength is plotted for a constant value of acceleration, then the 
maxima of the curves obtained for different values of acceleration 
lie at the Infinitesimal amplitude resonance frequency for high 
accelerations as mentioned above, and move to progressively 
lower frequencies as the piston acceleration is decreased. 

For the weak finite amplitude region where there are no 
shocks. Fig. 4-2 leads us to define resonance as that frequency 
for which the per cent of second harmonic is a maximum for 
constant piston acceleration. With this definition the tendency 
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for the resonant frequency to become lower as the acceleration is 
decreased continues to follow the same pattern as in the weak 
shock region. Unfortunately data were not collected to indicate 
whether or not this definition is consistent with the maximization 
of the average rectified voltage. The general trend for the 
resonant frequency to move to lower frequencies appears to 
continue as the infintesimal amplitude region is apprt>ached, 
however, no data were taken to follow this resonant frequency 
drift into the infinitesimal amplitude region. 
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CONCLUSIONS 



The most significant result of this study is the establishment 
of the ratio of acceleration to attenuation constant as a similarity 
parameter for waveforms in tubes of different cross-sectional 
area and shape. Once this ratio is known in two dissimliar tubes, 
their waveforms will be found to be the same for any value of 
acceleration as long as this ratio and the phase angle off 
resonance are the same. This similarity parameter is predicted 
theoretically for infinitesimal amplitude standing waves, but it is 
shown that the same parameter is applicable for weak finite 
amplitude and shock waves. The similarity parameter was 
investigated for only two types of tubes, but there is no reason to 
expect different behavior for different shapes of tubes with 
radically different attenuations assuming that the systems are 
"tight" in all cases. 

In order to determine th(.s similarity parameter, five methods 
of measuring attenuation In the infinitesimal amplitude region have 
been developed and have, been demonstrated to be . consistent withi 
each other and with theory. The values of attenuation constants 
obtained agreed with the results of other experimenters. These 
methods can be effectively adapted to the measurement of 
attenuation in a variety of gases under a wide range of 
pressures and temperatures. 
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In the region of weak finite amplitude waves the resonant 
frequency may be defined as the frequency causing greatest 
harmonic distortion for a given acceleration. This resonant 
frequency is observed to be lower than the resonant frequency 
for infinitesimal amplitude waves. This difference between the 
two resonant frequencies decreases as acceleration is increased. 
These observations are consistent with the results in the shock 
region where resonance is defined as the frequency for maximum 
shock strength. As the acceleration is increased in the shock 
region, the maximum shock strength occurs closer and closer to 
the resonant frequency of infinitesimal amplitude waves. 
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Figure 1 .1 EQUIPMENT SETUP 



I 





Figure 3-1. Variation of pressure magnitude at the rigid end as frequency is varied, maintaining tube 
length fixed. 
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varied, with frequency fixed at 200 cps. 



PHASE ANGLE IN DEGREES 




Figure 3-3. Variation of phase angle 0 as frequency is varied by 
an amount . ^f about the center frequency fQ. The 
smooth curve represents theoretical variation of 0 with At. 
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Figure 3-4. Comparison of attenuation measured by A/nP method 
with theoretical value in cylindrical tube . 
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Figure 3-7. Asymptotic approach of nPmax to a constant for 
high values of n. 
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^ Figure 3-8. Asymptotic approach of experimental value of 

I attenuation to theoretical value for high values of n. 
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PRESSURE MAGNITUDE 




Figure 3-9. Typical semi-log plot of exponential decay of pressure, 
Illustrating attenuation by reverberation technique. 
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Figure 3-10. Comparison of attenuation measured by metliod with 
theoretical value in rectangular tube. 
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Figure 4. 1 
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PEAK PISTON ACCELERATION IN METERS/SECOND 
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Figure 4-2. Second 'harmonic distortion of weak finite amplitude 
standing waves 
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Determination of Attenuation at 89 cps Resonant Frequency 





Cylindrical Tube 
(2.22 cm radius) 




1 

Rectangular Tube 

(3.05 X 6.90 cm 
Cross sectio 


n) 


Method 


Attenuation 


the or. 


Attenuation 


<A. theor. 


Theory 


1.17 X lO""^ cm"^ 




1.23 X lo"'* cm“^ 




Phase Angle 


1.34 .20 


1.15 


1.56 + .20 


1.26 


Bandwidth 


1.28 + .15 


1.09 


1.37 + .15 


1.11 


Reverberation 


1.29 + .10 


1.10 


1.45 + .10 


1.18 


Ratio of 
^min ^max 


1.76 + .50 


1.51 






Ratio of 
Acceleration 

^max 


1.17 + .20 


1.00 


1.37 + .20 


1.11 



Table 1 

(See Appendix II for error analysis) 
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APPENDIX I 



THEORY OF STANDING WAVES IN A CLOSED CYUNDER DRIVEN 

SINUSOIDALLY 



The values of attenuation measured in the round tube were 
found to be highly dependent on the tightness of the system. An 
excellent way to minimize leakage and at the same time to 
minimize alignment problems would be to seal both ends of the 
tube rigidly and to shake not Just one piston but the whole system. 
In this case the only possible source of leaks would be a small 
gasket where the probe microphone would be inserted. 

Assuming again a dissipative plane wave and applying the 
boundary conditions a = at x = 0 and x = L we 

obtain 

A - ^s) 



from which 



B -- 
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Applying these results to the expression for pressure previously 
obtained in section (3. l|and evaluating at x = 0 and x = L 
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Rewriting ^ 



in terms of transcendental functions 




n 



I - (^e.oSh »<L Sl^h e<t 



S/h h»{l c*s KL CosK S)'n HL 



K [" I . ^sin^,*kL Cts Hi^^^cosh^k Sm 

Again making the assumptions previously made. 



P 



p r r (hcos »Ly Y*<t. 

L <^0J i- s^h^^rl J 

This expression approaches a maximum for values of KL = (2n-l)TT 
and a minimum for KL = 2n TT 
Evaluating P at KL = (2n-l)TT 



p =. a/^ 

<^K L 

= 2. A 

pc?n- 0 fr 

This expression should provide an experimental method of 
obtaining ©< for n different values of KL. The coefficient 
may also be found by the bandwidth method. The same result, 
o( = ^ K holds for this system as for the system with only 
one end driven. 

The advantage of having a very tight system, characteristic 
of this system, was offset by several disadvantages. First, the 
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system was weight limited, which limited the size of tube to be 
shaken and eliminated the possibility of water jacketing the tube. 
Second, the length of the tube could not be varied. Third, and 
most serious , the microphone mounting could not be satisfactorily 
isolated from the vibration of the system. The microphone was 
extremely sensitive to vibration and the sensitivity was a very 
critical and unpredictable function of frequency. 

The prediction that pressure would be at a maximum at odd 
integral values of n was verified. The prediction that pressure 
would be at a minimum for even integral values of n was not 
verified, since relatively weak but still discemable resonance 
peaks were found there. These weak peaks are probably the 
result of harmonics in the driver input and do not invalidate the 
theory. The only value of c< which was closely consistent 
with theory was that obtained from the reverberation technique, 
which indicated that a tight system had been obtained, but not 
significantly tighter than the system with only the driver piston 
moving. Under the circumstances, further investigations of this 
technique iiid not seem warranted. 
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APPENDIX II 



ERROR ANALYSIS FOR ATTENUATION MEASUREMENTS 

1. Phase Angle Method. Phase angle measurements may be in 
error by several degrees due to drift in the vertical position of the 
oscilloscope trace. Frequency measurements may be in error by a 
maximum of 0 . 2 cps . Although determination of from a single 
measurement of 0 and ^f may vary by as much as 50% from the 
average value, averaging six readings equally spaced about the 
center frequency should reduce the error to less than 20%. 

2. Ratio of Acceleration to Pressure Method. This method is 
dependent on the accuracy of accelerometer and microphone. Assum- 
ing a calibration uncertainty of 1 db in the microphone and 0.5 db in 
the accelerometer, the overall systematic error is estimated at 15%. 

3. Bandwidth Method. As in the phase angle method the major 
source of error is in the determination of frequency, which should 
not be in error by more than 0.2 cps. Since the bandwidth was 
about 1.5 cps maximum error should be approximately 15% for 
frequencies around 100 cps. 

4. Ratio of P , to P Method. System noise due to 60 cps hum 

min max 

increased the value of P^in estimated 20%. Although this 

method is dependent on a ratio of pressures which should make it 
independent of calibration errors, the large difference between the 
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two pressures makes it necessary to read the corresponding voltages 
on different scales of the VTVM, thus introducing a possible error in 
voltmeter calibration. 

5. Reverberation Method. Since the driver piston is stopped for 
this method the boundary condition for this end is not the same as 
for the other methods, so the actual attenuation is not necessarily 
the same. The spread of data points on a semi-log graph for a 
single determination of c\ is approximately 3% but separate 
determinations of vary as much as 10%, apparently due to 
changing end conditions. 

An additional problem in comparing these methods for 
consistency is that it is extremely difficult to duplicate end 
conditions for a series of attenuation determination, even with a 
well-designed rigid-end piston. A certain determination of o(^ may 
be entirely valid, but if the microphone piston is withdrawn from 
the tube and re-inserted the attenuation may change due to the 
piston being in a slightly different position. 

It should be kept in mind that the method of determining c<, 
which gives results closest to the theoretical value is not 
necessarily the best method. Investigations of attenuation have 
generally indicated values of attenuation 12 to 18% higher than the 
theoretical prediction. 
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